
*** Verzeichnis festlegen
cd "C:\Users\...\replication_for upload\analysis for 1996-2021\replication_for upload\analysis for 1996-2021"

use allbus96-21.dta, clear

tab rd01
rename rd01 konfession
tab konfession, nol
tab konfession

*** konfessionslose in 4 codieren 
*recode konfession (6=4) (4=6)

label list rd01

*label define rd01 4 "KONFESSIONSLOS", modify
*label define rd01 6 "AND CHR RELIGION", modify


* kirchgangshäufigkeit

tab rp01
tab rp01, nol

gen kirchgang_häufig = 0
replace kirchgang_häufig = 1 if rp01 < 4


tab pt06
rename pt06 vertrauen_rK
rename pt07 vertrauen_eK
rename pt05 vertrauen_bundeswehr
rename pt08 vertrauen_justiz
rename pt03 vertrauen_bundestag

tab pa04
tab pa04, nol

tab pv01
tab pv01, nol
gen union_wahl =0
replace union_wahl = 1 if pv01==1

gen afd_wahl =0
replace afd_wahl = 1 if pv01==42

gen gru_wahl =0
replace gru_wahl = 1 if pv01==4

gen spd_wahl =0
replace spd_wahl = 1 if pv01==2

gen fdp_wahl =0
replace fdp_wahl = 1 if pv01==3

gen union_pi = 0
replace union_pi = 1 if pa04==1


gen rk_kirchgang = 0
replace rk_kirchgang = 1 if konfession==3 & kirchgang_häufig==1

gen rk_no_kirchgang = 0
replace rk_no_kirchgang = 1 if konfession==3 & kirchgang_häufig==0


gen ev_kirchgang = 0
replace ev_kirchgang = 1 if konfession==1 & kirchgang_häufig==1

gen ev_no_kirchgang = 0
replace ev_no_kirchgang = 1 if konfession==1 & kirchgang_häufig==0

gen fev_kirchgang = 0
replace fev_kirchgang = 1 if konfession==2 & kirchgang_häufig==1

gen fev_no_kirchgang = 0
replace fev_no_kirchgang = 1 if konfession==2 & kirchgang_häufig==0

gen oB_kirchgang = 0
replace oB_kirchgang = 1 if konfession==6 & kirchgang_häufig==1

gen oB_no_kirchgang = 0
replace oB_no_kirchgang = 1 if konfession==6 & kirchgang_häufig==0


* Geschlecht
tab sex
tab sex, nol
gen female = 0
replace female = 1 if sex==2 

* Alter

tab age

* 

gen konfession_kirchgang = 0
replace konfession_kirchgang = 1 if rk_kirchgang==1
replace konfession_kirchgang = 2 if rk_no_kirchgang==1
replace konfession_kirchgang = 3 if ev_kirchgang==1
replace konfession_kirchgang = 4 if ev_no_kirchgang==1
replace konfession_kirchgang = 5 if fev_kirchgang==1
replace konfession_kirchgang = 6 if fev_no_kirchgang==1
replace konfession_kirchgang = 7 if oB_kirchgang==1
replace konfession_kirchgang = 8 if oB_no_kirchgang==1

tab konfession_kirchgang


exit
* analyse konfession

* Modell 1
eststo clear
eststo: logit union_wahl i.konfession##year ///
female age
*margins, at(year=(1980 1984 1988 1992 1996 2000 2004 2008 2012 2016 2018 2021) konfession=(1 3 6))
*marginsplot, level(95) x(year) recast(line) xlabel(1980 1984 1988 1992 1996 2000 2004 2008 2012 2016 2018 2021) 
margins, at(year=(1996 2000 2004 2008 2012 2016 2018 2021) konfession=(1 3 6))
marginsplot, level(95) x(year) recast(line) xlabel(1996 2000 2004 2008 2012 2016 2018 2021) 

* Modell 2; mit Kirchgangshäufigkeit


* Union, KIrchang*Konfession, ab 1996
* replace konfession_kirchgang = 1 : RK + kirchgang 
* replace konfession_kirchgang = 2 : Rk + no kirchgang
* replace konfession_kirchgang = 3 : EV + Kirchgang
* replace konfession_kirchgang = 4 : EV + no kirchgang


eststo: logit union_wahl i.konfession_kirchgang##year ///
female age

margins, at(year=(1996 2000 2004 2008 2012 2016 2018 2021) konfession_kirchgang=(1 2 3 4))
marginsplot, level(95) x(year) recast(line) xlabel(1996 2000 2004 2008 2012 2016 2018 2021) 


esttab using "C:\Users\...\model_union.rtf" , replace ///
	abs b(3) se(3) ///
	star(+ 0.10 * 0.05 ** 0.01) ///
	aic ///
	scalars("ll Log pseudolikelihood") ///
	nonumbers /// 
	coeflabels(_cons "Constant" ///
	female "Female" ///
	age "Age") ///
	nonotes addnotes("Note: Numbers in parentheses are standard errors. Significance levels: + = p < 0.1; * = p < 0.05; ** = p < 0.01.") 

	